Single Server Systems — II. Busy Periods 

By S. 0. RICE 

(Manuscript received July 10, 1961) 

This is the second of two papers dealing with single server systems. 
Statistical problems associated with the busy periods, i.e., the periods during 
which the server is continuously busy, are considered in the present paper. 
The input and the service time distributions may be quite general, and the 
length of the waiting line is unrestricted. Among the new results is an 
asymptotic expression for the probability density of the lengths of the busy 
periods. This expression holds when the arrival rate is almost equal to the 
service rate and when the busy periods tend to be long. It is hoped that the 
methods used here will throw additional light on known results. 

I. INTRODUCTION 

This is the second of two papers dealing with single server systems. 
In the first paper, two subjects were discussed, (i) a method which 
led to the average values of several quantities of interest and (ii) the 
statistical behavior of a single server loss system. The present paper is 
concerned with the busy periods (the periods during which the server is 
continuously busy) in a single server system when no restrictions are 
placed on the queue length. 

The distribution of the busy period lengths has been studied by a 
number of investigators, among them E. Borel, 2 D. G. Kendall, 3 F. 
Pollaczek, 4 ' 5 L. Takacs, 6 ' 7 and B. W. Connolly. 8 ' 9 

A closely related problem is that of the storage of water behind a dam. 
An interesting survey of this subject has been given by J. Gani. 10 The 
moving server problem treated by McMillan and Riordan 11 and by Kar- 
lin, Miller and Prabhu 12 is also related to the busy period problem. Again, 
for Poisson input, the distribution of the busy period lengths is related to 
the distribution of the delays in "last come, first served" type of service 
(see Riordan and the references to earlier work given there). 

The most general results are those due to Pollaczek. His work leads to 
the joint distribution function of n,S,y where n is the number of services 
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comprising the busy period, S the length of the period, and y the length 
of the following idle period. 

The object of the present paper is to obtain, in another way, results 
equivalent to those derived by Pollaczek. Integral equations similar to 
Lindley's 13 equation for the waiting time distribution are first set up. 
These equations are then solved by methods similar to the Wiener-Hopf 
technique used by W. L. Smith 14 to solve Lindley's equation. The aim of 
the presentation is to illustrate the methods; the discussion is heuristic; 
and no claims are made for completeness or rigor. On the other hand, it 
is hoped that the different point of view will add something to the 
understanding of the earlier results. 

In Section II a brief review of some of the earlier results is given. In 
Section III the length of the busy period, as measured by the number n 
of customers served, is considered. The results of Section III are gen- 
eralized in Section IV to apply to the joint distribution of n, S, and y. 
Special attention is paid to the distribution of S. The transient behavior 
of the queue length is considered in Section V. Since the method used 
here to study the busy times is closely related to Smith's method of 
dealing with waiting times, a review of the waiting time problem is 
given in the Appendix. Some changes are made in Smith's method in 
order to make it fit in better with the busy problem. 

Throughout the paper, except possibly for (2), the intervals between 
arrivals are considered to be independent of each other and of the serv- 
ice times (recurrent input). The service times are also taken to be inde- 
pendent of each other and of the input. 

As mentioned in the companion paper, I am indebted to John Riordan 
for much help in the preparation of this paper. I am also indebted to V. 
Bene§ and L. Takacs for helpful discussions and a number of references. 

II. PRELIMINARY REMARKS 

In this paper, A(t) will denote the distribution function for the inter- 
vals between successive arrivals, B(t) the service time distribution, and 
A'(t),B'{t) the respective probability densities. The arrival rate is a, 
the service rate is 6, the average interval between arrivals is 

a' 1 = f [1 -A(t)]dt 

and the average service time is 

b' 1 = f [1 - B(t)]dt. 
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The Laplace-Stieltjes transforms defined by 
a ( s ) = / e ~ st dA(t) 

p(s) = f e~"dB(t) 

Jo 



(1) 



play an important role. 

Some information on the average length of a busy period (chosen at 
random from the universe of busy periods) may be obtained from the 
general results of the companion paper. 1 Thus, from Table I of that 
paper, when the service rate b exceeds the arrival rate a so that statisti- 
cal equilibrium exists, the average length S of a busy period and the 
average number n of customers served (during a busy period) are given 
by 

o _ ar 1 _ I , . 

b~^a~bW(0J' n ~W(6)' {2) 

Here the system capacity N is infinite, t is the average length of an idle 
period, and W(0) is the chance that a customer is served immediately 
upon arrival (zero waiting time). For Poisson input t is 1/a and 

5 * - 1 a /o\ 

s = b=- a - »=— ,• p = s (3) 

irrespective of the service time distribution. For inputs other than 
Poisson, one must know either t or W(0) in order to compute S from 
(2). 

For the sake of orientation we state the following known results. 

i. For Poisson input, general service and a < b, Kendall 3 has shown 
that the transform 

y (z) = f e - !S dG(S) (4) 

Jo 

of the distribution function G(S) of the busy period lengths satisfies the 
functional equation 

y(z) - 0[z + a - ay(z)}. (5) 

Here /3(s) is defined by (1). For Poisson input and exponential service, 
this leads to the probability density 

G'(S) = d ^l = b(SV^br 1 h(2SVab)e- (b+a)s (6) 
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where h(z) denotes a Bessel function of the first kind for imaginary 
argument. 

ii. Suppose that we are given a busy period which has just begun. 
Let/„ be the chance that it will end after exactly ft services. For Poisson 
input, Takacs 6 has shown that the generating function 

/(*)-Z>7. < 7 ) 

i 

must satisfy the functional equation 

f(x) = xfta - af(x)]. (8) 

For exponential service this leads to 

(2n - 2)! r l (9) 

Jn n! (ft- 1)!(1 + P) 2 "- 1 KJ 

and for constant service time (Borel ) to 

/.- W 1 ^- do) 

When a < b, all busy periods end eventually and 2"/» = *• When 
a > b, the customers arrive faster than the server can handle them, and 
sooner or later a busy period will start and never end. Given a busy 
period which has just begun, the probability that it will never end is 
1 — y}i°/« , where now E"/» * s ^ ess tnan 1 - For p °i sson input, expo- 
nential service, and a > b, the probability that the busy period will not 
end is 1 — (b/a). 

Hi. As mentioned in the Introduction, the results most closely related 
to the work of the present paper are those due to Pollaczek. 4 ' 5 In particu- 
lar he has shown the following. Given a busy period which has just 
begun, let G»'{y,S)dy dS be the chance that it will consist of exactly ft 
services, have a length between S and 8 + dS, and be followed by an 
idle period whose length lies between y and y + dy. Let 

y n (s,z) = ("dS rdye° v -* s G n '(y,S) (11) 

''O •'o 

where e v appears instead of e~ sv for later convenience. Then, for rather 
general input and service distributions, 

E x y n (s,z) = 1 - exp |- ^ ^ j (12) 

where 0gi|l, Re(«) < Re(f ) < < Re(z), Re(2: + f ) > 0, and 
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the path of integration for f runs from — ice to -H 00 in the strip speci- 
fied by the foregoing inequalities. In writing (11) and (12) it has been 
convenient to change from Pollaczek's notation to a notation resembling 
that used by Smith. 14 

III. NUMBER SERVED IN A BUSY PERIOD 

3.1 Derivation of Integral Equation 

Consider the busy period to start with the arrival of customer num- 
number 1 at time 0, and let the service time of the rth arrival be s r and 
t T the interval between the arrivals of customers r and r + 1. The busy 
period consists of one service if Si < ti , i.e., if iii = «i — t\ < 0. It con- 
sists of two services if 8\ ^ h and Si + So < h -\- t 2 , i.e., if U\ ^ and 
U\ + v-> < 0. In general, n demands (customers) are served in a busy 
period if Ui , U 2 , • • • , t/„_i ^ and U„ < where 

U n — V>i + U% + • • • + U n 

(13) 

U r = S r — t r ■ 

Since the w r 's are independent random variables whose distributions 
are known (indeed, avc cxp(— su r ) = /3(s)a(— s)) one may, in prin- 
ciple, find the joint distribution of Ui , U* , • • • , U n . Then the proba- 
bility/,, that n demands are served may be obtained by integration over 
the proper region in U\ , ■ ■ • , U n space. However, it is more convenient 
to use an indirect method which depends upon the solution of an integral 
equation similar to that for the waiting time distribution. 

Let pi(V)dV be the probability that V < lh < V + dV, and p n ( V)dV 
the probability that Ui , • • • , U„-i ^ and V < U n < V + dV, where 
V may be either positive or negative. Then 

Pl (V) = C'(Y) 

r , (14) 

Pn+i(V) = / Pn (v)C'(V -v) - 



dv 



where C'(t) = dC(t)/dt is the probability density of u r = s r — t r . 
Equations (14) lead to 

./(:»•,!') = xC'(V) + x [ J{x,v)C'(V - v) dv (15) 

J(x,V) = Zx n Pn (V) (16) 
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where the series and the integral are assumed to converge for ^ x £ 1. 
Equation (15) is the integral equation which must be solved to obtain 

J(x,V). 

The integral of p n ( V) from V = — « to gives the probability /„ that 
n demands are served and hence 

fix) = I>7 n = f J(x,V)dV. (17) 

When a < b, we expect all busy periods to end. Therefore 

th = l= f J(hV)dV. (18) 

1 •*— 00 

On the other hand, when the arrival rate a exceeds the service rate b 
we expect the queue length to increase with time, on the average. At 
the beginning of operations there may be a few busy periods, but even- 
tually a busy period starts which does not end. This state of affairs is 

indicated by the fact that S?/„ is less than * when a > b - Tnus » S iven 
a busy period which is just beginning, the probability that it will never 

end is 



i -Z/n = i - f J(i,v)dv. do; 

1 J— 00 



The length of the idle period following a busy period of length n 
(services) is — U„ . The chance that a busy period will be of length n 
and will be followed by an idle period whose length lies between y and 
;// + dy is p„(-y)dy. This is true for all values of a/b. When a < b, the 
chance that the length of an idle period picked at random from the 
universe of idle periods lies between y and y + dy is 

Zpn(-y)dy = J(i,-y)dy. (20) 

1 

3.2 Solution of Integral Equation 

The procedure used in the Appendix to solve Lindley's integral equa- 
tion (95) suggests writing J(x,V) as the sum of J-{x,V) and J+(x,V) 
where J_ is for V ^ and J+ is for V < 0. Multiplying (15) by 
exp( — sV) and integrating with respect to V from -co to + oo gives 

$_(.r,.s) +$ + (z,s) = a#(s)a(-*)[1 +*+(*>«)] ( 21 ) 

where Re(s) is confined to some suitable range and 4>+ is the Laplace 
transform of J + and <i>_ is similarly related to J_ [cf. (97)1. 
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Assume functions ty+(x,s), ^_(.r,s) and positive numbers Di , D* 

(which may depend on x) may be found such that 

.^( S W- s )-l=|±g£> (22) 

where 

(i) ¥+(x,8) is analytic in s and free from zeros in Re(s) > Di ^ 0, 
(ii) ^-(x,s) is analytic and free from zeros in Re(s) < D 2 with D 2 > 

JDi when < x ^ 1 and 
(Hi) 

ty+(x,s) —> s as [ s | — >■ oo in Re(s) > Di 

^_(.r,s) — > — s as | s | — > oo in Re(s) < Do. 

It is seen that & + (x,s) and ^_(.r,s) reduce to ^+(s) and \f/-(s) of the 
Appendix when x = 1. 

When .T/3(s)a(— s) is eliminated from (21) and (22) one obtains 

$_¥_ - ^_ - s = *+¥+ + ^+ - s (23) 

where s has been subtracted from both sides in order to keep them finite 
at infinity. Considerations of analytic continuation and the analogy 
with (101) suggest that both sides are equal to some quantity K(x) 
independent of s. When this turns out to be the case 

= gfr) - **M + s] 

V+(x,s) 

Suppose for the moment that 6 > a. To determine K( x) note (1) 
that <J> + (;c,s) is analytic in Re(s) ^ 0, and (2) that ty+(x,$) has a zero 
at s = .so = s (.r) where s a = when x = 1 and so w (1 — x)/(aT 1 — 
b~ ) when x is near 1. Also note that s is positive when b > a and x is 
slightly less than 1. Statement (1) is true since $ + (.r,s) is the Laplace 
transform of J+(x,V) and the integral 

/ J+{x,V) (IV = x Pr (fA ^ 0) + a: 2 Pr (17, , Us ^ 0) 

+ ••• < x(l - x)- 1 

converges. Statement (2) follows from relation (22) and the fact that 
the expansion of xP(s)a(—s) — 1 in powers of s shows that it is zero 
at s= so . Hence either ty + (x,So) is zero or^_(.r,s ) is infinite. The second 
possibility is ruled out if 1 — x is small enough to make s < Do since 
^_(.r,.s) is analytic in Re(s) < Do. Thus K(x) must be taken to be 
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—so in order to keep *+(.r,s) from having a pole in the region Re(s) ^ 
at s = so . 

Equating the left-hand side of (23) to -s and solving for <*>_ gives 

£ <TV_U-,F) dV = •_(»,) = 1 + £=^ (24) 

J_(^F) = 1 p '" e '^_(^) ds (25) 

47T4 Jc— too 

where c is chosen so that the path of integration in the s-plane passes 
to the left of the singularities of <f>_(z,s). Comparison of (24) and (17) 
leads to 

E*V. = *-W»= 1-^-1+^. (26) 

In this result, <&-(x,s) and ¥+{x,s) are obtained from (22); and s , a 
function of x, is that zero of x/3(s)a(-«) - 1 which approaches s = 
as a; — > 1 . In the foregoing definition of s , b is supposed to exceed a. 
When a and b are in any ratio, the statement is amended to read, "s = 
So is the one and only zero of <tr+(x,s) in Re s > when ^ X < 1." 
When a ^ 6, so tends to as x — > 1, and when a > b, s tends to a posi- 
tive number as x — » 1. These statements about s are made on the 
strength of the examples given below in Section 3.3 and hence cannot 
be regarded as proved in general. 

When a < b, (26) gives $_( 1,0) = 1 in agreement with (18). When 
a > b, so is not zero for x = 1, *_(1,0) is less than 1, and (19) shows 
that the probability that a busy period will not end is 

1 - *- (1 > 0) - dm (a > b) (27) 

The relations corresponding to (24) and (25) for J+(v,V) are 

j\->rj +{x ,V) dV = * + (x,s) = -I + ^ (28) 

j + ( X]V ) = -1. f" '°° e ,r *+(a*0 ds (29) 

I'M •'c-ioo 

where c is chosen so that the path of integration passes to the right of 
the singularities of <£ + (.t,s). 



SINGLE SERVER SYSTEMS — II 



287 



3.3 Examples 

The special cases used as examples in the Appendix to illustrate the 
waiting time distribution will be used here for the busy period length. 
It has been pointed out earlier that ty±(x,s) reduces to \p±(s) for x = 1. 

Example a. Poisson Input, Exponential Service. Equation (22) to de- 
termine the ty's becomes 



-& + (x,s) _ xab — (a — s)(b + s) 



where 



(a- 


»)(& + «) 








s 2 + (b - a)s - (1 - 


v)ab 


(s — s )(s — 
(a- s)(6 + 


Si) 




(a- s)(b + s) 




■) 


3 


= a -^ ± W(b 


+ a) 2 


— 4akc. 





(30) 



(31) 



As a; runs from to 1 the roots s and Si move along the real axis in 
the s-plane as shown in Fig. 1 . 
The conditions for (22) lead us to take, as in (109), 



*+(x,s) = 



(s - gp) (a - Si) 
(« + b) 



^_(.r,s) = a — s 



(32) 



with Di — s and Z) 2 = a. Expressions (26) give the generating function 
for the probability /„ that, given a busy period which has just begun, 
it will consist of exactly n services. Setting s = in (32) to obtain 
SP-(x,0) leads to 



L *7- = i - 



So 



¥_(*,0) 



= 1-* 



a 
a + b 



(33) 



~ ?r V(b + a) 2 - 4ate. 
2a 2a 



The coefficient of x n in the power series expansion of the last expression 
a<b a>b 



1 = 

-b 



a-b 





Fig. 1 — Runges for s and s t as x runs from to 1. 
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gives the value (9) for /„ . When x = 1, the sum £"/" is equal to 1 
when a < b and to b/a when a > b. Thus, from (19) or (27) the proba- 
bility that the busy period will not end is 1 — (b/a), with b < a. 
Inserting the expression 

*_(*,«) = 1 + S -^° = ^^-° (34) 

a — s a — s 

obtained from (24) in the integral (25), and evaluating the integral 
by closing the path of integration on the right when V < 0, shows that 

{(a - s )e aV (7<0) 



J-(*,V) = 



(35) 
(V^O) 



From (20), the probability density for the lengths of the idle periods 
when a < b is 

J.(l-y) =ae- a!/ 

as expected for Poisson arrivals. When a > b, the idle period probability 
density is still given by the expression on the right, but now J-(\, — y) 
has to be divided by the normalizing factor 1 — (s /a) in which x = 1, 
i.e., by b/a, the probability that the busy period will end. 

Although $+(x,s), J + (x,V) are not needed to calculate the probabili- 
ties given above, their values as obtained from (28) and (29) will be 
stated for the sake of completeness 

s — Si 

f(6 + Sl )e'* v (V ^ 0) 
J+(x,V) = 

JO (7 < 0) 

Note that b + s x = a — s Q . 

Example b. Poisson Input, General Service (Takacs 6 ). In this case, 
(22) and the related conditions are satisfied by the analogue of (110) 

V+(x,s) = s - a + axp(s), ty-(x,s) = a - s (36) 

with D\ = so and D 2 = a. With the help of Rouche's theorem it may be 
shown that ¥+(x,s) has only one zero, s = s , in Re s > when ^ 
x < 1. Furthermore, since ^+(x,0) is negative and ty + (x,a) is nonnega- 
tive, < so ^ a. Equations (24) and (25) show that $-{x,s) and 
./_(.t,F) are given by (34) and (35), just as for Poisson input and ex- 
ponential service, and we still have 

y x »f n = i _ i° . (37) 

a 
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However, So no longer has the simple form (31), but it still tends to 
if a ^ b and to a positive number if a > b as x — » 1. 

When Lagrange's expansion theorem is applied to 1 — (s /a) it is 
found that (see Pollaczek, 5 p. 102, Eq. (8.37)) 



/„ - (_a) 



n\ 



[£" w >*]~ (38 > 



which gives (9) and (10) as special cases. It should be recalled that from 
(3) the average number served in a busy period, for Poisson input and 
a > b, is n = 1/(1 — p). 

Example c. Recurrent In-pat, Exponential Service. When one takes the 
steps leading to (111) (which pertains to the waiting time distribution 
for this case) as a guide, he is led to 

*+(M) = ( * " ?}* ~ *> ■ »..(*,)- |'r»" , 7* ) (39) 
b -\- s xba{—s) — b — s 

where Si is the only zero of 

h(s) = s + b - xba(-s) (0 < x < 1) (40) 

which lies in Re(s) > and s is the left-most zero of h(s) in Re(s) > 
when x is close to 1 . The existence of Si may be established with the help 
of Rouche's theorem. Then h(0) = (1 - x)b > and h(-b) = -xa(b) 
> show that — b > si > 0. To make the existence of s plausible, 
consider the case when a is nearly equal to b and x is close to 1 . When s 
is small, the series (114) for a(—s) gives 

h(s) = (1 - x)b + s(l - xba~ l ) - ^|?i + • • • ■ (41) 

It is seen that h(s) has a double zero at s = when x = 1 and a = b. 
When x = I — t, with e small and positive, and 1 — baT 1 = -q is small, 
the double zero splits into two simple zeros given approximately by 



= -eb - s V + 



ba 2 s" 



The two roots of this equation are small, real, and of opposite sign. The 
positive root corresponds to s and the negative one to Si . 
From (26) 

± x »f n = i + O-ag* = , + <}jl*>k (42) 

i ^+(.r,0) si 
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and Lagrange's expansion theorem gives (see Pollaczek, 5 p. 103, Eq. 

(8.42)) 

/» = C„_i - C n (Co = 1) 



C„ = 



ST£'^<-^>1 



When a < b, the average number served in a busy period may be ob- 
tained either by differentiating (42) and setting x = 1, or using the value 
of T^(0) obtained just below (111). Both methods give 



n = 



W(0) 



r--i • 

L sijx-i 



(44) 



Exam-pie d. Erlangian Input, Erlangian Service. Let f(s) denote the 
same polynomial in s as in Example d of the Appendix. Then 

•*•«-> - 1 - o + &r i 1 - sF - i - =# • 

When x = the polynomial a: — /(s) has a zero of order /c at s = — bk 
and a zero of order I at s = a£. Now let < rt < 1. From Rouche's 
theorem and the fact that \f(s) | > .r on the imaginary s axis, it may 
be shown that x — f{s) has k zeros in Re(s) < and I zeros in Re(s) 
> 0. Denote the zeros in Re(s) < by Si , • • • , s k , the left-most zero 
(when x is close to 1) in Re(s) > by s , and the remaining zeros in 
Re(s) > by s k+l , • • • , s k +i-i . Then (22) takes the form 

x - f(s) _ (s - so) (s — Si) • • • (s — Sfc+t-i) _ ¥+(a,a) ,^ . 



f(s) (s - al) l (s + bk) k *-(**) 

and the conditions on the ^'s are satisfied by 

_ , s (s — s )(s — Si) • • • (s — s k ) 

* +M = — ™\ — (46) 

, T , / v (s - aZ) 

^-(.TjS) = —7 r 7 r • 

(S — S fc+ i) • • • (S — S k+ i-.i) 

From (26) the generating function for the number served in a busy 
period is 

T* »* i ( n~> i (1 -x)(-bk) k 

2_,X f„ = 1 — {al) 8o8 k+ i8k+t • ■ • Sk+l-i = 1 — 

I Si • • • Sk 

where s , Si , • • • , s k+ i-i are functions of .r. 
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IV. DURATION OF BUSY PERIOD 

The distribution of the busy period lengths for recurrent input and 
general service may be obtained by an extension of the foregoing analy- 
sis. A number of the steps will be omitted here because of the similarity 
with the work of Section III. For the sake of simplicity, throughout this 
section it will be assumed that a < b and hence that all busy periods 
eventually end. 

4.1 Derivation of Integral Equation 

Let p n (V,S)dV dS be the probability that U 1 , • • • , U n -i ^ 0, V < 
U n < V + dV and S < S n < S + dS where U„ is defined by (13) 
and S n is the sum Si + s 2 + • • ■ + s n of the first n service times. It may 
be shown that p„(V,S) is zero for V > S and 

Pl (V,S) - B'(S)A'(S- V) 

p, l+l (V,S) = [ da- J dvp n (v,a)B'(S - a)A'(S - V - a + v). 
•'o *'o 

These may be combined to give the integral equation 

J(x,V,S) = xB'(S)A'(S - V) 

r s f , , (47) 

+ * / da / dvJ(x,V,<r)B(S - a)A(S - V - <r + v) 
Jo J u 

where 

J(x,V,S) = £,x n p n (V,S). (48) 

n=l 

The probability that a busy period will consist of exactly n services, 
have a length between »S and S + dS, and be followed by an idle period 
of length between y and y + dy is 

p n (-y,S)dS dy = G n '{y,S)dS dy (49) 

where G„'(y,S) is the density function introduced in Section II. The 
probability (hat a busy period will consist of exactly n services and 
have a length between S and »S + dS is dS times 

On(S) = f Pn (V,S)dV. (50) 

J— 00 
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Summing with respect to n shows that the probability of a busy period 
length between S and S + dS is dS times 

G'(S) = f J(l,V,S)dV. (51) 

J— 00 

Furthermore, the probability the busy period consists of exactly n serv- 



ices is 



/,,= ('dS f P n(V,S)dV. 

Jn -I— m 



4.2 Solution of Integral Equations 

Multiplying (47) by exp[— zS - sV], and integrating V from - qo 
to S and 8 from to » gives 

*_(x,s,2) + $+{x,sjs) = xj9(s + «)a(-s)[l + 3> + (.t,s,z)] 

*+( w) - j[" dS £ rfF e-* s -V(a;,y,S) (g2) 

<S>_Cw) = I" dS f dVe- a -' y J(x,V,S), 

From the last equation 

I e-' r J(x,V,S) dV = J-. I e !S *-(x,s,z) dz (53) 

where c is such that the singularities of the integrand lie to the left of 
the path of integration. 
Assume the factorization 

^« + ««-)- 1 -^ <54) 

where the ^ functions satisfy conditions (t), (»), (*«*) set forth in con- 
nection with (22), 2 being regarded as an imaginary constant (or is at 
least on the path of integration Re(z) = c in (53)). When z = 0, 
* ± (.t,s,z) reduces to the ^ ± (.r,s) of (22) and when z = and .r = 1, 
to the il/±(s) of (99). The analogue of (24) turns out to be 

, / v ,1 s- s (x,z) ,__>. 

*-(x,s,z) 
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where s = s (x,z) is that root of 

x${s + z)a(-s) -1=0 (56) 

which tends to s = as x — > 1 and z — » 0. 

The various probabilities of interest may be computed, at least in 
theory, from $-(x,s,z). The relations are indicated in Table I. The 
Laplace transforms y n (s,z), y n (z), y(z) are the same as those introduced 
in Section II. 

From (55) and the last two entries in the table it follows that "the 
Laplace transforms of 2 x"G n '(S) and G'(S) are respectively 



<M.i-,0,z) = 1 - 



so(x,z) 
*-(*Az) 



= 1 + 



T(.) -•-(«*)- 1 - i^ - 1 + 



[1 - xp(z)]s (x,z) 
* + (z,0, 2 ) 

[1 - (3(z)]s (l,z) 



* + (l,0,z) 



(57) 
(58) 



Furthermore, the generating function (26) for /„ may also be written 
as M z,0,0). 

4.3 Examples 

The special cases used earlier will again serve as examples. The results 
of Examples b and c are equivalent to results given by Pollaczek. 



Table I 



Probability Function 


Laplace Transform 


G n '(y,S) = Pn (-y,S) 


y n (s,z) 


= Expression (11) 

= f dS f dVe-**-'Vp n (V,S) 

•>0 J— to 


£ x»G n '(y,S) = J(x,-y,S) 
l 




00 

^2 x n y,,(s,z) = *_(x,s,z) 

l 


G n '(S) 


I? 


-' s G n '(S) dS = y n (0,z) = y n (z) 


£ x"G n '(S) = / J(x,V,S) dV 

1 J-00 




n 

£ x"y„(z) = +-(x,0,2) 
1 


G'(S) = / J(1,V.S) dV 

J 00 




-* S G'(S) dS = 7 (z) = *-(l,0,z) 
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Example a. Poisson Input, Exponential Service. In this case (54) be- 
comes 



^+(x,s,z) _ s 2 + s(z + b — a) — a(z + b — xb) 
ty-(x,s,z) (s + b + z)(a — s) 

(s — s )(s — Si) 



(59) 



(s + & + z)(a - s) 
where s> = Sj(x,z) is given by 



± i V(b + a + z) 2 - ixab. (60) 



So\ _ a — b — z 

Sij 2 "~" 2 

When x is fixed and z runs from — ico to -H 00 , So traverses an oval- 
shaped path in the s-plane. The oval lies in Re(s) ^ 0, it starts and 
ends at s = a (corresponding to z = ±t'oo ), and when z = it crosses 
the real s-axis between s = and s = a as indicated by Fig. 1, case 
a < b. At the same time, Si traverses a path roughly parallel to the 
imaginary axis. The path lies in the region Re(s) < 0, is asymptotic to 
the line Re(s) = — b at z = ±i«>, and crosses the real s-axis between 
s= —b and s = a — b. As x — > 0, the oval shrinks to the point s = a 
and the path of Si tends to the straight line Re(s) = —b. 

This behavior of So and Si and the similarity with (32) lead us to take 

* + (w) = (s ~ So)( _f ~ Sl) , *_( W ) = a - t (61) 
s -\- o -]- z 

From (57) 

#_(z,0,z) = 1 - ^^ = ^ [& + a + 2 - \/(& + a + z) 2 --4za6]. (62) 
a Za 

Inverting this Laplace transform (see, for instance, Pair 556.1, Camp- 
bell and Foster 16 ) gives (Takacs 6 ) 

Z x n G n '(S) = S' 1 a/- «" (frfs)fl /i(S\/S5) (63) 

i K a 

where Ji(z) is a Bessel function of the first kind for imaginary argument. 
Putting x = 1 in (63) gives Kendall's 3 expression (6) for G'(S). 

Example b. Poisson Input, General Service (Prabhu 17 ). In analogy 
with (36) 

ty+(x,s,z) = s — a + xa0(s + z), &-(x,s,z) = a — s (64) 
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and (57) gives 

$-(x,0,z) = 1 - a _1 so0,2) (65) 

where s = s (x,z) is that root of 

s - a + xa@(s + z) = (66) 

which tends to zero as x — > 1 and z — > 0. Rouche's theorem may be used 
to show that if ^ # < 1 and Re z ^ 0, then (66) has one and only 
one root (namely, s a (x,z)) in Re(z + s) > 0. 

When s in (66) is replaced by s (x,z) = a — a<f>_ (from (65) with 
$_ = $-(x,0,z)) one obtains 

$_ = x${z + a - o*_). (67) 

This is a functional equation satisfied by the generating function 

•--£«>.(«). (68) 

Setting x = 1 gives Kendall's equation (5) 3 

7(2) = /3(2 + a - a 7 (z)) 

for the Laplace transform of G'(S). 

Example c. Recurrent Input, Exponential Service (Connolly, 8 Takacs'). 
The analogue of (39) is 

, v (s - s )(s - ft) 

* +W; S + 6 + 2 ' 

(69) 

*_< w) = ( ; zj^g - g) 

.roa(— 5) —8 — — 2 

where s = Si is the only root of 

xba(-s) - s — 6 — 2 = (70) 

which has a negative real part when Re(z) ^ and ^ x Si 1, and 
s = So is the root which tends to s = when 2 — > and x —* 1. 
Equations (57) and (58) give 

*_ ( „o )2) = 1 + « + (i-*)fc 

Sl (x,z) (n) 

7(2) = 1 + -7^-r 

which may be used in Table I. The coefficients of —2 and z'/2 in the 
power series expansion of 7(2) give the first and second moments of the 



296 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1962 

distribution of the busy period length S. With the help of (70) it is 
found that 

SlO . 4 

^_____ C7Q\ 
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where 5(0 is a unit impulse. Since the right-hand side of (75) is a con- 
tinuous function of S, it cannot hold in this case. However, formal cal- 
culation of (75) gives 

G'(S) m 8 - 3/2 (2Tbr !i e- bsll - p)212 . (77) 

By using n\ ~ \J2tu n n e~ n and p « 1 it may be shown that (77) be- 
comes a smoothed version of (76) when S is replaced by rib~\ A better 
approximation to (76) may be obtained by noting that the Laplace 
transform 7(2) of (76) is of the form ^/ n exp(— znb~ x ). Hence 7(2) 
is a periodic function of z with period *2t6, and contributions to (73) 
may be expected not only from the region around z = but also from 
regions near z = i2irkb where k = • ■ • , — 1, 0, 1, 2, • • • . It may be 
shown that the sum of these contributions does indeed approximate 
(76) when n is large and p ~ 1. 

The approximation (75) holds for large values of 8. One may obtain 
an idea of the behavior of G'(S) for small values of £ by noting that 
now the busy period contains only a small number of services. In particu- 
lar, the assumption that the very short busy periods consist of a single 
service leads to the approximation G'(S) &Gi'(S) = B'(S)[l - A(S)]. 
Higher-order approximations may be obtained by computing pi(V,S), 
P2(V,S), • • • step by step and using (50). The first step gives the result 
cited, namely 

Gi(S) = f B'(S)A'(S - V) dV = B'(S)[1 - A(S)}. (78) 
J— 00 



V. GROWTH OF QUEUE 

In Section III the functions ^ ± (x,s) were introduced to obtain the 
chance /„ that a busy period will consist of n services. It is interesting 
to note that these functions may also be used to determine W r (t), the 
waiting time distribution for the rth customer. Only a sketch of the 
procedure is given here. More complete information on W r {t) is given 
by Pollaczek, 5 Spitzer, 18 and Lindley. 13 
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so that L+(x,t) = for t < 0, it is found from (94) that 

L+(x,t) = WM + x £ L+(x,i - r) d('(r) . (80) 

This holds only for / ^ 0. To solve the integral equation, let L-(x,t) 
denote the value of the right-hand side for t < and set L-(x,t) = 
for t ^ 0. It is found that 

*-(*,«) + <P+(X,8) = s" 1 + x<p+(x,s)P(s)a(-s) (81) 

where <p + is the Laplace transform of L + and ^>_ is similarly related to 
L_ (cf. (97)). Equation (81) is similar to (21). When the functions 
ty±(x,s) appearing in (22) are introduced, the analogue of (23) turns 
out to be 

*_(*,«)*-(*,«) - ■ 4 *-(<M) + s"^_(x,0) 

= e> + (x,s)¥ + (a;,s) + s _l *_(*>0) 

where s _1 ^_(a;,0) is added to both sides in order to cancel the pole of 
s _1 ^_(.t,s) at s = 0. Setting both sides of (82) equal to a quantity K(x) 
independent of s leads to 

. [s - s (x)]9jxfl) 

where s = s (x) is the one and only zero of ^ + (x,s) in Re(s) > 0. 
Inversion gives the required result: 

± x -w r (t) = ^i i r +io ° [s - »/*y <** (e > o) . (83) 

V s (z) 2in J.-** s^ + (a;,s) 

5.2 Poisson Input and Exponential Service 

For the special case of Poisson input and exponential service, the val- 



ues of M> ± (z,s) are given by (32). Insertion in (83) gives 

Zx'Writ) -i / 

27TZ •'e-i 



1 f +,w q(a + 6)e" rfg 

S ° S(S " *> " (84) 

1 (a + b)e ' , ' /- T _ T v/i- CI 

(1 — VI — cx)e 



1 - x 26(1 - x) 
where e > 0, t > and 

4ab /(6 + a) 



■ ■ w+w - - — (85) 
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Upon using 



i - vr^) t+1 = E 



(cx) n+, (2n - k)\(k + 1)2* _1 



7=k (n + ])!(n - /c)!4" 
(1 - Vn^i-)exp[r-rVn^] - ~ £ T^TT^n P M 

ii n =0 (,»* "T 1 ^ II 

it is found that 

PT«) _ 1 ae- bt g CP(n,r) 
k M) a + bh?o{n + 1)14" 

where P(n,r) is the polynomial 



[86) 



(87) 



p , v ^ (2n-fc)!(fc + l)(2r) fc 

t=o /c!(w-fc)! ( 88 ) 

= (2r) n+5 eV- i [(l + r)/C 1+J (r) - r/C„_j(r)] 

and K denotes a Bessel function of the second kind for imaginary argu- 
ment. 19 It may be shown that P(n,r) is 0(4"w!/\/n) as n — > w with t 
finite. 

Setting x = 1 in (86) and using the result to transform (87) leads to 

W r (t) = 1 - p e- ( "- a)t 4- a <r"F r (0 (a ^ 6) 

TT r (0 = + a e- bl F r (t) (a ^ b) 

where 

The value of F,(<) is unchanged when the values of a and b are inter- 
changed. 

When a = b = 1, differentiation of the generating function (84) leads 
to expressions for the probability density W T '{t), r > 0: 

JF/(0 - 2^j (*) UC +5 (/) - /C r _j(0] 

. _e^_ y^ 1 (2r-n- 1)1(2/.)" 
2 2 --~o r\n\(r- n - \)\ K ' 



2 /"" cos Ittdttl" 1 £ ~| 

" jr A (w 2 + l) r [u 2 + 1 " 2rJ" 
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5.3 Queue Behavior as r — > » 

Return now to the case of recurrent input and general service. When 
r becomes large W,(t) approaches the limit W(t) discussed in the Ap- 
pendix if a < b. Questions related to this approach have been studied 
by Lindley, Pollaezek, Spitzer and others. 

When a > b, the customers arrive faster than the server can handle 
them and the waiting line tends to grow steadily. An idea of the behav- 
ior of W r (t) in this case may be obtained by using the notation of the 
Appendix. The rth customer arrives at time T T = to + k + ■ • • + t r -\ 
and his service begins at »S* r + I T , where S T = s + Sj -f- • • • + s r _i 
and I r is the total amount of time the server is idle in the interval (0,T r ). 
The waiting time of the rth customer is w r = S r + I r — T r . Since the 
server may be expected to be continuously busy after a few initial idle 
periods, we expect 7 r to approach some constant value as r — > °o . Thus, 
I r becomes small in comparison with S r — T r = u Q + U\ + • • • + u r -i ■ 
Upon making the approximation w T « S r — T r and using the central 
limit theorem, it is found that, for a > b and r » 1 



W r (t) = W — t= exp -—— . 

at a\/'2irr L *ro- J 



(90) 



Here /x = b~~ l — a -1 is the average value of u r = s r — t r and a' = v u + 
v A is its variance, v B and v A being the respective variances of s r and i r . 
When a = b, the approximation (90) no longer holds. However, 
Pollaezek 5 has shown that, at least for Poisson input, in place of (90) 
we have 



«>0). (91) 



W T '{t) ^—4-expl-^ 
This agrees with the asymptotic form of the integral in (89). 
Appendix 

Waiting Time Distribution* 

In Sections III and IV the busy period problem has been investigated 
by a method similar to the Wiener-Hopf technique used by Smith 14 to 
deal with the waiting time distribution. The application of the method 
to waiting time problems is reviewed in this Appendix. Smith's approach 
has been changed slightly in order to make it fit in better with the busy 
period problem. 

* Here the "waiting time" of a customer is the interval between his arrival and 
the instant his service begins. 



302 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1982 

A.i Derivation of Integral Equation 

The basic integral equation is due to Lindley. 1 Suppose the single 
server system (with unrestricted queue length, recurrent input, and 
general service) starts operations when the 0th customer arrives at time 
0. Following Lindley, write w r for the waiting time of the rth customer, 
s T for his service time, and t r for the interval between the arrivals of 
customers r and r + 1. The rth customer stays in the system for an 
interval of length w r + s r , and the (r + l)th customer arrives t r units 
of time after this interval begins. If w r + s r ^ t r , the (r + l)th cus- 
tomer does not have to wait and w T +i is zero. When w r + s r > t T , the 
(r + l)th customer has to wait w T + s r — t r units of time. Thus 

W r+ X = (W r + S r ^ t T ) 

(92) 

W r +1 = W r + S T - t r (w r + S r > t T ) . 

It is assumed that the s r are independent random variables with the 
common distribution function B(t), and also that the t T are independent 
(of each other and of the s r ) with common distribution function A(t). 
The rule for combining probability distributions shows that the distri- 
bution function C(l) of the variable u r = s r — t r is given by 

C(t) = / B(t -\- t) dA(r) . (93) 

Since w r and u T are independent random variables, it follows upon re- 
writing (92) as 

w r+ \ = (w T + u r ^ 0) 

Wr+l = W r + U r (w r + U r > 0) 

that the distribution function for w T+ \ is 

i(0 = f_ Writ - r)dC(r) (^0). (94) 



ir 



r+H 



By starting with 1F (0 = 1 for t ^ 0, which states that the 0th cus- 
tomer is served immediately, one may compute Wi(t), W 2 (t), • • • , in 
succession from (94). When the service rate b exceeds the arrival rate 
a, i.e., when p = a/b < 1, W T (t) tends to W(t) as r — > oo where W(t) 
satisfies Lindley 's integral equation 

W(t) = f W(t - r) dC(r) (t^0). (95) 

J- 00 

17(0 is the distribution function for the waiting time when statistical 
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equilibrium prevails. Note that (95) holds only for nonnegative /. When 
t is negative, W(t) is zero but the integral on the right does not vanish. 

A.2 Solution of Integral Equation 

Let W-(t) be the value of the integral for t < and take W..{t) = 
for t ^ 0. Since W(t — r) is zero for t > t, (95) may be written as 

WM) + W(t) = f W(t- t)C'(t) dr (96) 

which holds for all real values of t. The derivative C'{t) = dC(r)/dT 
is the probability density of the random variable u r = s r — t r . 

Multiply both sides of (96) by exp(-s/), where < Re(s) < D 
with D such that the following integrals converge. The existence of D 
is ensured by assumptions made below. Integrate from t = — <» to 
t = co and introduce the transforms 

*>+(«) = C c-'W(t)dt= f e~"W(t)dt 

J_oo Jo 

?_(s) - [ e" l W-(t) dt = ( c- l WM) dt (97) 

J— co J— co 

[ e~"C'{r) dr = ave exp [-ss r + st T ] = /3(s)a(-s) 

J— 00 

where /3(s) and a(s) are the respective Laplace-Stieltjes transforms of 
the service and interarrival distribution functions B{t) and A(t). This 
carries (96) into 

*>_(s) +<p+(s) = <p + (s)P{s)a(-s) 

*»_(«) = V+(s)\P(sM-8) - 1]. 

Since W(t) and B(t) are distribution functions, both <p+(s) and /3(s) 
are analytic in the region Re(s) > 0. To ensure convergence of the inte- 
grals in (97) involving W-(t) and C'(r), assume that the probability 
density A'(t) = dA(t)/dt exists and is 0[exp(— Dt)] as t — » <», where 
D is positive but may be arbitrarily small. It may then be shown from 
(93) that C(t) is 0[exp(D/)] as t — > - co and, using this in (95), that 
W-(t) is also 0[exp(Z>0] as t — > — «». It follows that both <p-(s) and 
a(—s) = ave exp(s^) are analytic in the region Re(s) < D. 

Now suppose that functions \p+(s) and ^_(s) may be found such that 

0(s)a(-s) -l-*4^ (99) 
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where (i) yp+(s) is analytic and free from zeros in the half-plane Re(s) > 
0, and (ii) ^_(s) is analytic and free from zeros in Re(«) < D. Although 
these functions may be expressed as integrals when suitable conditions 
are satisfied (see, for example, Smith 14 ) their expression in tractable 
form is usually the most difficult step in obtaining W{t). For future 
convenience assume that \p+(s) and i^_(s) may be chosen so that, 

\b+(s) — > s as I s I — > co in Re(s) > 

(100) 
\f/-(s) — -> — s as | s | — > co in Re(s) < D. 

The difference in sign is required by the fact that the left-hand side of 
(99) tends to — 1 as s — > ±z°° (unless both A(t) and B(t) have dis- 
continuous jumps, a case we shall rule out). 

When the resolution (99) is possible, (98) becomes 

*>_(*)*-(*) = *+(«)*+(«>, < Re(«) < D. (101) 

The right-hand side is analytic for Re(s) > and the left hand for Re(s) 
< D. Equality in the strip implies that each is the analytic continuation 
of a function which has no singularities in the finite part of the s-plane. 
It turns out that when conditions (100) are satisfied, this function may 
be taken to be a constant K. Indeed, conditions (100) were imposed to 
make this so. Then the Laplace transform of W(t) is 

K ' 

<p+(s) = -^- (102) 

which is analytic in Re(s) > by virtue of the requirements on \p+(s). 
Since 

lim s<p + (s) = lira f e~" dW(t) = I 

it follows that ^+(0) is and K is given by 

K = \unt±^l = \*t&>] = ^ + '(0) . (103) 

»-o s {_ as J.=o 

The constant K is also equal to W(0), the probability that a customer 
will not have to wait for service. One way to see this is to note that from 
(102) and (100) 



K = <p + (s)\f/+(s) = lim <p+(s)s = lim s I e "Wit) dt 

= lim / e-"W f- J du = W(0) . 



(104) 
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The limit exists since < W(t) ^ 1 and W(t) decreases monotonically 
as / decreases to 0. Although W(t) is discontinuous on the left at t = 
and may have discontinuities for positive values of t (as it docs when 
both .4(0 and B(t) have discontinuities), we take it to be continuous 
on the right at t = 0. 

Thus when ^+(s) is known, <p+(s) is determined and W(t) may be 
obtained by inversion, 

Wit) =£-.r m rr\ (oo). (105) 

The Laplaee-Sticltjes transform of W(t) is 

w ( s ) = / e-'dW(t) = $v+(s) = -r-y-v. (100) 

As s — » oo , w(s) tends to IF(0) = K. The value of -dw(s)/ds at s = 
is cental to the average waiting time. 

In some cases ^_(s) is simpler than tf+(s) and in place of (102) one 
may use 

K_ 

[0(s)a(-s) - !]*_(«) 



*+(«) = r^.x / .7 rrrvrr ( 107 ) 



where differentiation of (99) gives 

K = * + '(0) = (a -1 - 6 _1 )^_(0). (108) 

a. 3 Examples 

The results just obtained will be illustrated by several special cases, 
all of which have appeared in the literature. 

Example a. Poisson Input, Exponential Service. Here a(s) = a/ (a + 
s), 0(s) = 6/(6 + s), and ^+(s), £-(«) are to be determined from (99), 
which now takes the form 

ba _ s(s + b — a) _ ip+(s) 



(b + s)(a - s) (6 + s)(o-s) *_(*)' 

Inspection shows that when D = a 

f + (,) = s(s +_^~ fl) , *_(«) = a - s (109) 

satisfy the requirements set forth in connection with (99). Then (103), 
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(b -a) 

= 1 - p 



(102) and (105) 


give 




K = 


-- lim 

«-0 


s 



<P+(s) = 



b 
K (1 - p)(& + s) 1 



^+(s) s(b — a + s) s b — a + s 

w(t) = i - P <r (6 - a) '. 

Incidentally, for this case the probability density of u r = s r — t T is 

[ab(a + fc)"V" (m < 0) 

C"(«) = 

[a5(a + b) _I e" 6u (w ^ 0) 

Example b. Poisson Input, General Service. In this case (99) becomes 

POQa _ , = s — a + qff(s) = ^+(s) 
a — s a — s ^-(s) ' 

With the help of Rouche's theorem (Titchmarsh, p. 116) and the fact 
that [1 — /3(s)]/s is the Laplace transform of 1 — B(t), it may be shown 
that s — a + a/3(s) has no zero in Re(s) > when b > a. 
Then, with D = a, 

\ft+(s) = s - a + afi(s), ^-(s) = a - s, 

K = (a -1 - 6 _1 )^_(0) = 1 - p = W(0), (110) 

stf s(l - p) 



w(s) = s^-+(s) = 



^+(s) [s - a + a/3(s)] ' 



The expression for w(s) is sometimes called the Pollaczek-Khinchin 
formula. The coefficient of — s in the power series expansion of w(s) is 
equal to the average waiting time. 

Example c. Recurrent Input, Exponential Service. Replacing /3(s) by 
b/(b + s) in (99) gives 



ba(—s) — b — s \~ba(— s) — b — s \ \s(s — Si 



b + s ! s(s — Si) J [_ b + S 



D 






where Si is the one and only zero of s + b — ba(—s) which lies in Re(s) 
< when b > a. The existence of s% may be established with the help 
of Rouche's theorem. 1 If s = is the only new zero of s + b — ba(—s) 
introduced when the region Re(s) < is extended to Re(s) < D, one 
may take 

( v s(s - si) , v s(s - Si) ( v 

^ + (s) = —— , ^_( s ) = -— . (Ill) 

b + s ba{—s) — b — s 
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These functions satisfy the conditions stated just below (99). Then 

K = lim Mil = * = W(0) 
«-o s b 

( \ K s.(6 + s) 1 _ 1 + si&~' 

v+{8) ^ + (s) == bs(s - si) s s - Si 

W(t) = 1 - (1 +s 1 6"V 1 '. 

Example d. Erlangian Input, Erlangian Service. In this case the input 
and service time probability densities are 

,/ /lN al(alt) _1 -„(/ „'/.>. bk(bkt) -bkt 

A M = irnTi e ' B {l) = w^w e 

and have the Laplace transforms 

It is found that (99) becomes 

/(a) M«) ' JW V at) \ bk) 

For a < b (the only case considered here), the polynomial F(s) = 
1 — f(s) has a zero at the origin, zeros Si , ■ ■ • , s* in Re(s) < and 
zeros Sfc+i , • • • , s*+i_i in Re(s) > 0. This may be shown with the help 
of Rouche's theorem. It turns out that si , • • • , s k lie inside a circle of 
radius bk centered on s = -bk, and s k +i , • • • , s k+ i-\ He inside a circle 
of radius al centered on s = al. Hence 

. , v s(s — Si) • • • (s - s t ) 

*+<•) = — (s + bkr — ' 

*_(,) = - ( « ~ ° ' r (112) 

(S - Sa+i) • • • (S - S/H-l-l) 

(— bk) k s/t+i • • • s fc+ ,_i 

For the case of regular arrivals and constant service time, inspection 
shows that W(t) = 1. It appears that in this case letting k = 2 — > *> 
should give 
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Example e. Arrival Rate Almost as Large as Service Rate. It often 
happens that when the arrival rate a is almost as large as the service 
rate b, most of the customers have to wait a long time for service. In 
such cases one may obtain an approximation for the waiting time dis- 
tribution function W{t). 

Let Oi , a 2 , hi , &a be the first and second moments, and v A , v B be the 
variances (which are assumed to exist in the following discussion) asso- 
ciated with A(t) and B(t). Then a\ = 1/a is the average spacing be- 
tween arrivals and bi = 1/6 is the average service time. From the inte- 
gral (105) for W(t) it is seen that the behavior of W{t) for large values 
of t is determined by the right-most singularities of K/\f/ + (s). In Exam- 
ple a (Poisson input, exponential service), these are poles at the zeros 
of yp+(s) which occur at s = and s = Si = a — b. Note that «i is nega- 
tive and S\ — > as a tends to b. Furthermore, the value of ^-(s) near the 
origin does not change markedly as a — > b; and the same is true for the 
remaining factor (6 + s) _1 appearing in ^+(s). 

Many other cases show the same type of behavior as a — > b. In gen- 
eral, the function (99) 

has a double zero at s = when a = b. When a becomes slightly less 
than b, and a\ slightly greater than b\ , one of these zeros remains at 
s = and the other moves to s ~ s x where 

si = ■ < 0. (113) 

V A + Vg 

This may be seen upon using 



0(8) = 1 - h* + ^ + 0(S 2 ) 



a( —s) = 1 + «is + -^- + o(s 2 ) 



(114) 



and assuming that (a t — W) 2 is negligible in comparison with v A + v B . 
The approximation (113) for Si is given by Smith who points out its 
importance in the present case. 
An examination of the earlier examples leads us to take 

4,+(s) « s{s - Si)C (115) 

when s is near the origin. Here C is a constant equal to the value of 
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the remaining portion of ^+(s) at s = 0. Equations (103) and (105) 
then give 

K = hmt+^1 = - Sl c (116) 



K -s t 1 



1M«) s(s - si) s s - si (117) 

If (0 « 1 - e' 1 '. 

Thus, when a —> b and v A + v B does not tend to zero, W(t) is given by 
the approximation (117) where Si is given by (113). The average wait- 
ing time is — \/s\ . It should be noted that this approximation gives 
W(0) « instead of the true (small) value W(0) = K = -SiC. Un- 
fortunately, there appears to be no simple expression for C correspond- 
ing to (113) for s t . 

a. 4 Conclusion 

Finally, it will be mentioned that when a customer departs after be- 
ing served, the chance p n that n customers remain in the system is equal 
to the chance that an arriving customer will find n in the system. Fur- 
thermore, for ^ x < 1 



2- •<• Pn = 9 ■ / 

h=0 ATI J ,— i 



a(-s)p(s) Kds 



-,» [I - .ra(-s)] \J/ + (s) 

where c > and is such that the singularities of /3(s)/^ + (s) and a( — s)/ 
[1 — xa{— s)] lie on opposite sides of the path of integration (this re- 
stricts A (t) and B{t) somewhat). The right-hand side may be replaced 
by one of several integrals which differ slightly from the one shown. 
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